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ABSTRACT 

‘‘This  paper  develops  three  distinct  models  for  studying  perishable 
inventory  systems.  The  perishable  items  have  a deterministic  usable  life 
after  which  they  must  be  outdated.  For  each  of  the  models,  analytical 
expressions  are  found  for  steady-state  distributions  which  characterize 
the  inventory  systems.  Knowledge  of  this  steady-state  behavior  may  be 
used  for  evaluation  of  system  performance,  and  for  consideration  of 
alternatives  for  improving  system  performance.  — 

S^The  first  model  considered  assumes  that  both  the  demand  process 
and  the  inventory  replenishment  process  are  stochastic  processes  that 
may  be  modelled  as  Poisson  processes.  The  second  and  third  model  assume 
that  inventory  is  replenished  by  a constant  production  process.  The 
second  model,  assuming  continuous  inventory  units,  has  Poisson  demand 
requests  with  the  size  of  each  request  distributed  as  an  exponential 
random  variable.  The  third  model  has  Poisson  demand  requests  with  all 


demands  being  for  a single  unit 
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1.  INTRODUCTION 

This  paper  is  concerned  with  the  development  of  models  for  studying 
inventory  systems  in  which  the  item  is  perishable.  In  particular,  the 
item  is  assumed  to  have  a lifetime  of  m periods,  after  which  it  is  said 
to  be  outdated  and  cannot  be  issued  to  service  demand.  Examples  of  such 
items  might  be  blood  in  a blood-banking  system,  packaged  food  items  or 

chemical  supplies.  This  paper  considers  three  distinct  models  of  perish- 

% 

able  inventory  systems.  For  each  of  these  models,  a steady-state  distri- 
bution of  the  inMeutory  status  is  found.  Knowledge  of  this  steady-state 
behavior  may  be  used  for  evaluating  a system's  performance,  and  for 
considering  alternative  means  for  improving  the  system. 

Much  of  the  literature  studying  perishable  inventory  systems  has 
dealt  with  finding  optimal  periodic  ordering  policies.  Typical 
of  this  work  is  the  work  by  Nahmias  [10],  [11],  [12],  Fries  [4],  and 
Cohen  [5].  Implicit  in  this  work  is  the  assumption  that  an  order  can 
always  be  filled  by  the  supplier.  The  analysis  necessary  for  finding 
optimal  policies  requires  an  m-diraension  state  space  where  m periods 
is  the  lifetime  of  the  item.  Consequently  this  work  has  dealt  with 
determining  and  evaluating  approximate  policies.  The  most  common  policy 
considered  is  an  order-up-to  policy:  each  period  an  order  is  placed  so 
to  bring  the  total  Inventory  up  to  a preset  level.  Related  to  this  work  is 
that  of  Jennings  [7]  who  used  a simulation  model  to  study  a blood  banking  system. 
With  the  simulation,  he  was  able  to  evaluate  various  ordering  policies, 
and  was  able  to  consider  the  effects  of  centralized  control  on  a regional 
multi-hospital  blood  banking  system. 

The  models  developed  in  this  paper  are  distinct  from  the  work  on 


periodic  ordering  policies.  In  that  here  specific  replenishment  policies 
are  assumed.  The  replenishment  or  ordering  policy  is  not  a decision 
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variable;  rather^  given  a replenishment  policy,  the  models  are  descriptive 
of  the  stochastic  system  behavior.  In  this  sense,  this  paper  is  related 
to  the  descriptive  Markov  models  for  perishable  inventory  systems  presented 
by  Pegels  and  Jelmert  [13],  Brodheim,  Derman  and  Prastacos  [1],  and  Chazan 
and  Gal  [2] • 

Pegels  and  Jelmert  [13]  model  a blood-bank  inventory  by  a Markov  chain. 

This  work  has  been  criticized  by  Kolesar  [8].  Brodheim,  et  al.  study 
an  inventory  system  assuming  a constant  replenishment  policy;  that  is, 
each  period  n units  of  new  inventory  are  supplied.  They  model  the  system 
as  a Markov  chain  in  discrete  time,  but  are  unable  to  solve  the-  steady- 
state  equations.  They  do  establish  upper  and  lower  bounds  on  various 
measures  of  interest,  such  as  expected  outdates  and  expected  shortages, 
provided  that  n,  the  supply  quantity,  is  less  than  the  expected  demand 

i 

ii 

rate.  Chazan  and  Gal  consider  an  inventory  system  operating  under  j 

! 

an  order-up-to  replenishment  policy.  They  use  a Markov  model  to  find  j 

1 

j 

bounds  on  the  expected  outdates  assuming  a general  demand  process.  For 
Poisson  demand,  they  are  able  to  find  the  steady-state  behavior  for  the 
system  for  the  Markov  chain  in  continuous  time.  That  is,  whenever  a 
unit  is  requested  or  a unit  becomes  outdated,  it  is  assumed  that  a new 
unit  is  instantly  acquired  to  replace  the  old  unit  to  maintain  a constant 
Inventory  level;  note  that  for  such  a system  there  can  never  be  any 
shortages.  Uslig  the  results  from  the  continuous  case,  tight  bounds 
are  developed  for  the  Markov  chain  in  discrete  time. 

The  remainder  of  the  paper  is  organized  as  follows:  the  next  section 
considers  an  inventory  system  for  which  the  replenishmen*’  process  is 
stochastic  and  may  be  modeled  as  a Poisson  process.  Sections  3 and  A 
consider  models  in  which  there  is  constant  replenishment.  The  problem 
in  these  sections  is  similar  to  that  considered  in  [1);  however,  here 
we  are  able  to  present  analytical  results.  Section  3 considers  a continuous 

L ' 
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model  in  which  the  size  of  demand  requests  is  exponential*  Section  4 
assumes  that  demand  requests  are  all  for  single  units*  In  all  models  it 
is  assumed  that  requests  arrive  as  a Poisson  process.  The  final 
section  gives  a summary  and  discussion  of  the  results  of  the  paper. 
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2 . STOCHASTIC  REPLENISHMENT 

Consider  the  inventory  system  assuming  that  both  the  demand  process 
and  the  replenishment  process  is  Poisson,  Demand  requests  occur  at 
rate  y,  and  a request  is  always  for  one  unit.  The  inventory  is  replenished 
at  rate  A with  unit  replenishments.  This  type  of  replenishment  process 
may  be  appropriate  whenever  there  is  uncertainty  in  the  supply  procedure. 
For  instance,  in  the  case  of  blood  inventory,  to  a certain  extent  supplies 
depend  on  volunteer  donors  who  may  be  thought  of  as  behaving  as  a random 
process.  Alternatively,  in  a production  environment,  the  production 
process  may  be  very  sensitive  to  human  or  random  elements,  and  consequently 
the  time  to  produce  a successful  unit  is  variable.  An  example  of  this 
might  be  a repair  depot  where  the  time  to  diagnose  the  failure  or  problem 
is  random,  and  consumes  most  of  the  repair  time.  In  addition  to  stochastic 
demand  and  replenishment,  assume  that  the  item  has  a lifetime  of  m periods, 
the  issuing  policy  is  to  issue  the  oldest  usable  item,  and  that  demand 
requests  that  cannot  be  filled  from  Inventory  are  turned  away  or  handled 
exogenous  to  the  system. 

This  system  can  be  modeled  as  a Markov  chain  in  continuous  time. 

One  characterization  of  the  state  space  is  by  a complete  account  of  the 

inventory.  That  is,  the  state  variable  would  be  {N(t),  A^(t) , 

where  N(t)  is  the  number  of  units  in  inventory  at  time  t,  and  is 

th 

the  age  of  the  i unit  in  inventory  where  the  units  are  ranked  by  age, 

for  i«l,2, . . . ,ll(t) . In  this  form,  however,  the  analysis  of  the  process 
is  difficult,  if  not  impossible,  due  to  the  varying  dimensions  of  the  state 
space. 

As  an  alternative,  consider  the  process  A(t)  where  A(t)  is  the  age 


of  the  oldest  unit  in  inventory.  Given  the  Poisson  assumptions  for  demand 
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and  supply,  and  assuming  we  have  rw  other  information  about  the  inventory 
status,  then  this  process  is  Markov  and  will  provide  a sufficient  charac- 
terization of  the  system  for  our  purposes.  The  process  A(t)  can  take  on 
any  value  a for  0 £ a < m.  The  state  space  for  A(t)  must  also  Include 
the  state  when  there  is  no  inventory  in  the  system;  this  will  be  denoted 

by  A(t)  = E.  Note  that  A(t)  = E is  distinct  from  A(t)  = 0. 

The  transitions  for  A(t)  can  be  characterized  by  examining  the 
interval  (t,  t+A)  for  the  following  three  cases,  where  A > 0 is  small. 

Case  1 - Suppose  that  A(t)  « a for  0 ^ a < m-A.  Then  at  t+A,  we  have 

that  either  a demand  request  occurred  (with  probability  yA)  or 
it  did  not  [with  probability  (l-yA) ] . If  no  demand  occurs, 
then  A(t+A)  = a+A.  If  a demand  occurs,  then  the  oldest  unit 
(age  = a)  is  issued,  and  A(t)  becomes  the  age  of  the  next  oldest 
unit  in  inventory  or  A(t)  = E if  the  issued  unit  was  the  only 

unit  in  inventory.  Define  x to  be  a random  variable  for  inter- 

arrival time  of  units  into  inventory;  T is  an  exponential 
random  variable  with  mean  1/A.  If  x < a+A,  then  there  is 
a second  unit  in  inventory  and  A (t+A)  = a+A-x.  If  T ^ a+A, 
this  implies  that  the  oldest  unit  was  the  only  unit  in  inventory, 
and  the  system  is  now  empty,  A(t+A)  = E.  Note  that  for  this 
case  a replenishment  epoch  has  no  effect  on  the  process. 

Case  2 - Suppose  that  A(t)  =*  a for  m-A  ^ a<m.  Then  at  t+A  the  oldest 
imit  will  have  been  issued  or  outdated.  Thus,  using  similar 
reasoning  as  in  Case  1,  A(t+A)  = m-x  for  X <_  m,  A(t+A)  = E 
for  X > m. 

Case  3 - Suppose  that  A(t)  = E.  In  (t,  t+A),  the  system  will  either 


remain  empty  [with  probability  (1-AA) , A(t+A)«E]  or  a new 


L 
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replenishment  will  arrive  [with  probability  AA,  A(t+A)=0]. 

Given  the  exponential  interarrival  times  for  replenishments, 
the  probability  of  a replenishment  in  (t,  t+A)  is  independent 
of  the  time  of  the  last  previous  replenishment.  Note  also  that 


this  case  is  unaffected  by  demand  requests,  since  they  are  not 

■k 


backordered. 


Considering  the  above  three  cases,  the  process  {A(t)}  is  clearly 
Markov  since  the  transition  laws  for  the  interval  (t,  t+A)  depend  only 
on  the  state  A(t)  at  time  t. 

Define  p(x,t)  to  be  the  probability  density  for  A(t)  = x,  0 £ x < m, 
and  let  TrCt)  be  the  probability  that  ACt)  = E.  The  equations  of  motion 
can  be  shown  to  be 


(1) 


3p(x.t)  3p(x.t) 

3x  dt 


= -yp(x,t)  + u p(y,t)Ae  + p(m,t)Ae 

y X 


0 < X < m. 


and 


(2)  ^.(t)  = 


-XirCt)  + p p(y,t)e  + p(in,t)( 


-Am 


y=0 


In  addition  to  this,  the  total  probability  mass  must  be  1: 


(3)  TT(t)  + p(x,t)dx  = 1. 


By  differentiating  both  sides  of  (3)  with  respect  to  t,  and  substituting 
(1)  and  (2),  we  find  a boundary  condition  for  the  process: 


* 

The  current  model  is  easily  extended  to  allow  backorders  provided 
that  A > p. 
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(4)  p(0,t)  = XTT(t). 

The  steady-state  equations  are  found  by  letting  t p(x,t)  -*•  p(x), 

and  7T(t)  -*■  7T.  Equations  (1),  (3),  and  (4)  are  now  rewritten  as 

(5)  ^ ^ -Vip(x)  + y p(y)^e  + p(ni)Ae 

0 ^ X < m, 

(6)  IT  + p(x)dx  = 1, 

and 

(7)  p(0)  = Xtt. 


The  solution  to  (5)  , (6)  , and  (7)  is  fo'vind  to  be 


(8)  p(x)  = Ke 


(X-u)x 


for  0 < X < m 


(9) 


7T  = K/X 


(10)  K = 


provided  that  X y.  For  the  special  case  when  X = n,  the  solution  is 

(8),  (9),  and  ; 

i 

j 

(11)  K = X/(mA  + 1).  i 

i 

i 

From  (8)  - (10)  we  can  now  express  analytically  four  measures  of  i 

system  performance  where  an  outdate  occurs  when  a unit  reaches  age  m and  ■ 

cannot  be  issued,  and  a shortage  occurs  when  a demand  request  is  turned  j 

away  because  there  is  no  inventory. 

a)  the  expected  number  of  outdates  per  time  unit  (0) . 

b)  the  expected  number  of  shortages  per  time  unit  (S) . 

iA 
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c)  the  expected  age  of  an  Issued  unit  (A). 

d)  the  expected  number  of  units  in  inventory  (I). 


Using 

simple  probabilistic  arguments,  we  have 

(12) 

0 = 

p(m) 

(13) 

S » 

yir 

(14) 

A = 

xp(x)dx/^^Q/“  p(x)dx 

(15) 

I = 

p(x)(l+Ax)dx 

The  expressions  (12)  - (15)  can  be  used  to  evaluate  a system,  given 
the  system  parameters  X,  y,  and  m.  These  expressions  can  also  be  used 
to  examine  the  sensitivity  of  the  measures  to  the  system  parameters. 

For  instance,  taking  A = 1*0,  y = 1.0,  m = 20,  as  a base  case 

Figures  1,  2,  and  3 show  the  effect  of  varying  the  replenishment  rate, 
the  demand  rate,  and  the  unit  lifetime,  respectively. 

The  present  model  assumes  that  demand  requests  arrive  as  a Poisson 
process,  and  are  for  single  units.  The  model  can  be  easily  extended 
to  the  case  where  the  size  of  the  demand  request  is  distributed  as  a 
geometric  random  variable  d.  That  is,  for  d being  the  number  of  units 
demanded , 

(16)  pr(d=n)  * (1-p)”  ^p  for  n=*l,2,... 

(17)  E(d)  » l/p 

The  analysis  of  this  problem  is  similar  to  the  previous  case  except  for 

differences  in  the  transition  laws  when  demand  requests  occur.  Define 
Tg  to  be  the  sum  of  d replenishment  Interarrlval  times,  where  each  inter- 
arrival  time  is  an  independent  exponential  random  variable  with  mean  l/X. 


SCALE  FOR  6,  S (ITEM  UNITS) 


FOR  A (TIME  UNITS) 
FOR  i UTEM  UNITS) 
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Since  d is  a geometric  random  variable,  x is  a compoimd  geometric  distri- 

s 

bution  of  exponential  random  variables.  It  can  be  shown  that  x^ 

is  an  exponential  random  variable  with  mean  1/pA.  Now  for  A(t)  « a < m - A, 

if  a demand  request  occurs  during  (t,  t+A) , then  A(t+A)  = a + A - x 

s 

if  X < a + A,  or  A(t+A)  = E if  x > a + A.  Note  that  x corresponds 
s s s 

to  the  removal  of  d units  from  invt-ntory.  However,  if  m-A  £ A(t)  < m, 
then  the  oldest  unit  is  outdated  during  (t,  t+A)  with  probability  1 - yA 
(i.e,  no  demand).  In  this  instance,  letting  X be  a single  replenishment 
interarrival  time,  we  have  A(t+A)  = m-x  for  X < m and  A(t+A)  = E for 
T > m. 

Using  the  same  reasoning  as  in  the  previous  case,  the  steady-state 
equations  for  this  system  are 


(18) 


dp(x)  ^ 
dx 


-yp(x)  + y „ p(y)  Xe  + p(m)  Xe 


y«x 


0 < X m, 


and  equations  (6)  and  (7),  where  X = pX.  The  solution  to  the  three 
equations  (6),  (7),  and  (18)  is 


(19)  p(x)  = (^)e-e 


(20)  7T 

(21)  K 


K 

X 


1 + 


(X-y)x  ^ |X-^  j ^Xx  ^-m(X+y-X) 

X-^\  -m(X+y-X) 


(¥) 


X-y-A  . (A-y)m  (X-t  . 1 


(^y  ■"A) 


-1 


Equations  (12),  (14),  and  (15)  may  be  used  to  determine  expected 
•k 

outdates,  expected  age  , and  expected  inventory  level.  The  expected 
shortages  per  time  unit,  however,  now  consist  of  shortages  which  occur 
when  there  is  no  inventory,  and  shortages  which  occur  when  a particular 


Here,  the  expected  age  needs  to  be  reinterpreted  to  be  the  expected 
age  of  the  oldest  unit  issued  to  fill  a demand  request. 
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deaiand  request  depletes  the  inventory  present*  The  component  of  shortages 
which  occur  when  the  system  is  empty  is  pTT/p,  since  the  expected  size  of 
a request  if  1/p  and  the  expected  number  of  requests  per  unit  time  when 
the  inventory  is  zero  is  yTT,  To  determine  the  expected  shortages  caused 
by  demand  requests  in  excess  of  present  (non-zero)  inventory,  suppose 
A(t)  = X and  a demand  request  occurs.  Letting  I be  the  total  Inventory 
level  at  time  t,  then  given  A(t)  = x we  have 


vn-1  -Xx 


for  n=l,2, . • , 


Combining  (22)  with  the  distribution  for  the  size  of  a demand  request 
(16),  the  expected  shortages  conditioned  on  A(t)  = x and  the  occurance 
of  a demand  request,  denoted  as  S(x),  is 


S(x)  = I I (ia-n)»pr[I=n|A(t)=x]*pr(d=ro) 
n=l  m=n 


Hence,  for  this  system,  the  expected  shortages  per  time  unit  may  be 
written  as 

(24)  S = M^/p  M p(x)S(x)dx 

o 

» y7T/p  + K[1  - e ^^](l-p)/p. 

where  the  first  component  is  shortages  when  the  system  Is  empty,  and  the 
second  component  is  shortages  from  partially-filled  orders. 


-14- 


3.  CONSTANT  REPLENISHMENT  WITH  EXPONENTIAL  DEMAND  REQUESTS 

Consider  an  Inventory  system  for  which  the  replenishment  process 
is  constant;  that  is,  inventory  is  continuously  produced  at  a constant 
rate  of  c units/ time  unit.  This  type  of  replenishment  is  found  in  any 
continuous  production  facility,  such  as  a chemical  plant  or  a refinery, 
where  it  is  very  costly  to  shutdown  or  interrupt  the  process.  For  this 
system,  it  is  assumed  that  demand  requests  arrive  as  a Poisson  process 
at  rate  \i  with  the  request  size  distributed  as  an  exponential  random 
variable  with  mean  l/y.  Implicit  in  this  model  is  the  assumption  that 
Inventory  may  be  treated  as  a continuous  entity.  In  addition,  it  is 
again  assumed  that  Inventory  is  issued  oldest  first,  inventory  expires 

it 

at  age  m,  and  inventory  shortages  are  not  backordered.  The  time  unit 
for  the  system  is  defined  such  that  the  production  rate  c = 1. 

Similar  to  the  previous  system,  this  inventory  system  may  be  charac- 
terized as  a Markov  process  with  state  variable  A(t)  corresponding  to 
the  age  of  the  oldest  unit  in  inventory.  Note,  however,  that  A(t)  now 
also  represents  the  amount  of  inventory  present.  The  process  {A(t)} 
ranges  over  the  interval  [o,m].  The  transitions  for  the  process  are 
described  by  considering  the  following  two  cases: 

Case  1 - Suppose  A(t)  = a,  0 £ a < m -A.  Then  A(t+A)  = a+A  if  no  demand 
requests  occur  in  (t,  t+A) . If  a demand  occurs,  then 
A(t+A)  « (a+A-x)^,  where  x is  the  exponential  random  variable 
for  the  size  of  the  demand  request. 

Case  2 - Suppose  A(t)  « a,  m-A  £.  ot  £ m.  At  t+A,  if  a demand  has  occured, 

7 


The  backorder  case  is  a simple  extension  of  the  current  model  provided 
that  y < Y* 
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then  A(t+A)  » Otherwise,  A(t+A)  « m.  Note  that  in 

r the  latter  case.  Inventory  is  being  continuously  outdated* 

f 

1 

i 

I Define  p(x,t)  to  be  the  probability  density  for  A(t)  “ x,  0 ^ x £ m, 

f and  7r(t)  to  be  the  probability  mass  that  A(t)  = m*  The  equations  of 

L 

I motion  for  this  system  are 

(25a)  |_  p(x,t)  + 1^  p(x,t) 

= -vip(x,t)  + yTT(t)  + /“  yp(y,t)Ye”^^^”’‘^dy, 

y*x 

for  0 ^ X £ m, 

and 

C25b)  = -yirCt)  +p(m,t). 

i In  addition,  we  have  the  following  boundary  conditions: 

I 

[ 

i 

I (26)  p(0,t)  « y-0'^”*  yp(y,t)e  ^^dy  + yTT(t)e 

1 

By  letting  t ->  »,  p(x,t)  p(x),  and  7r(t)  tt,  the  steady-state 

i 

\ equilibrium  equations  are  found  to  be 

i 

(27)  ^ p(*^  “ -yp(x)  + yiTYe'^^™”’'^  + yp(y)Ye"^^^~’‘^dy. 

for  0 £ X £ m, 

(28)  0 ■ -yiT  + p(m),  and 

(29)  p(0)  - y-o'^”  yp(y)e'^^dy  + yire”^. 

The  solution  to  (27)  - (29),  using  the  fact  that  total  probability 


mass  is  one,  is 


(30)  p(x)  = Ke 


(31)  TT  - -P(ib)!i 


-(u-Y)x 


0 < X < m, 


i^ere 


(32)  K 


provided  M ^ y.  For  M - Y.  the  solution  is  (30),  (31),  and 


(33)  K » y/d-hnp). 


By  probabilistic  reasoning,  we  can  state  analytic  expressions  for 
the  expected  outdates  per  time  unit  (0),  the  expected  shortages  per  time 
unit  (S) , the  expected  age  of  the  oldest  unit  supplied  for  a demand 
request  (A),  and  the  expected  inventory  level  (I): 


(34)  0 ® TT 


(35)  S - (p  p(0) 


(36)  A 


(37)  T 


xp(x)dx  + mir 


xp(x)dx  + mu 


Inventory  units  are  being  continuously  outdated  at  the  unit  rate  of 
production  (c*l)  whenever  A(t)  = m which  occurs  with  probability  mass  tt* 

To  find  the  expected  shortages,  if  a demand  request  is  larger  than  available 
inventory,  then  due  to  exponential  demand  sizes  the  expected  number  of 
unfilled  requests  is  l/y*  p(0)  reflects  the  frequency  of  occurance  of 
demand  requests  which  completely  deplete  the  inventory*  Note  also  that 
the  expected  age  of  the  oldest  unit  issued  equals  the  expected  inventory 
level;  this  is  true  since  the  age  of  the  oldest  unit  equals  the  Inventory 
level,  and  since  demands  occur  at  random  epochs*  Given  the  expressions 
(34)  - (37)  as  functions  of  the  system  parameters  p,  y»  model 

may  be  easily  used  for  system  evaluation  and  trade-off  studies. 
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4.  CONSTANT  REPLE2JISHMENT  WITH  UNIT  DEMAND  REQUESTS 

The  system  studied  In  this  section  Is  Identical  to  that  studied 
in  the  previous  section  except  that  all  demand  requests  are  for  one  unit 
of  inventory.  Again  the  system  may  be  characterized  by  the  process 
{A(t)}  for  A(t)  being  the  age  of  the  oldest  unit  in  inventory  at  time  t. 

To  analyze  this  system,  we  will  not  study  the  equilibrium  flow  equations 
for  the  system,  but  rather  will  draw  upon  results  from  queueing  theory. 

It  is  useful  to  think  of  the  system  as  a single-server  queue,  where  the 
production  process  is  the  server,  the  inventory  is  a positive  backlog  in 
the  queue,  and  a demand  request  is  a customer  arriving  at  the  queue 
requiring  one  unit  of  service.  In  this  framework,  the  system  may 
be  modeled  as  a capacitated  M|d|1  queue,  where  the  process  {m-A(t)} 
corresponds  to  the  virtual  waiting  time  for  the  m1d|1  queue.  When 
A(t)  » m,  the  server  is  idle  and  the  queue  empty;  for  A(t)  = ^ < m,  the 
server  is  busy  with  total  queue  length  Q equal  to  the  smallest  integer 
greater  than  m-A(t),  and  with  a virtual  waiting  time  of  m-A(t).  Since 
A(t)  is  defined  io  be  non-negative,  this  implies  that  the  queueing  system 
has  capacity  m;  queue  arrivals  which  would  put  the  virtual  waiting  time 
in  excess  of  m are  truncated  so  that  the  virtual  waiting  time  equals  m. 

Hence,  the  determination  of  the  distribution  for  {A(t)}  is  equivalent 
to  finding  the  distribution  of  the  virtual  waiting  time  for  a bounded 
M|d|1  system.  Unfortunately,  to  this  author’s  knowledge,  this  problem 
is  unsolved  (see  Gavish  and  Schweitzer  [5]  who  considered  a closely-related 
problem).  However,  results  exist  for  finding  the  distribution  for  a discrete 
approximation  of  {A(t)}.  The  approximation  is  analogous  to  using  queue 
length  as  a surrogate  for  the  virtual  waiting  time.  Here,  A(t)  = m 
corresponds  to  an  empty  queue,  ro-1  ^ A(t)  < m corresponds  to  a queue  length 


it 

This  is  appropriate  here  since  each  customer  requires  one  unit  of  service. 
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of  one,  and  so  on.  Defining  to  be  the  probability  that  n ^ A(t)  < n+1 
for  n**0,l,  • . • ,in-l  and  p^  to  be  the  probability  that  A(t)  “ m,  then  the 
results  of  Keilson  [8]  for  the  finite-capacity  m1g|1  system  can  be  extended 
to  the  current  problem  for  y < 1: 


(38)  p^  » 1 - (l-p^)/y 

(39)  p,  “ . i«l,2,...,m 

1 m— 1 

m 

(40)  K - [e^  + y(  Z e.)] 

i-0  ^ 

where  e^^  is  the  steady-state  probability  of  queue  length  i for  the  m|d|1 
queue  with  infinite  capacity.  That  is  (see  [6]) 


(41) 

(42) 

(43) 


1 - y 


(l-y) (e^-1) 


(1-y) 


1-1 


(-D^y^^e 


k ,k  y(i-k) 


k=0 


k! 


odsi:-  , 

k=0 


i-2  / ,sk„k  y(i-k-l) 


(i-k-l)‘ 


k! 


for  i=2,3, . . • 


Keilson  [8]  also  presents  a solution  for  y > 1;  however  this  result 
is  in  terms  of  a Green's  function  and  will  not  be  discussed  here. 

Now  (38)  - (40)  can  be  used  to  find  the  four  evaluative  measures  of 
system  performance: 


1 


(44)  0 - j 

1 

i 
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These  expressions  are  analogous  in  their  derivation  to  (34)  - (37), 
However,  the  above  measures  are  approximate  since  they  correspond  to  a 
discrete  representation  of  the  process  for  the  oldest  unit  in  inventory. 

The  system  considered  in  this  section  is  identical  to  that  considered 
by  Brodheim,  et  al.  [1];  however  the  results  are  quite  different.  In  [1] 
the  system  is  modeled  as  a Markov  chain  in  discrete  time  on  discrete 
state  space;  given  this  representation,  the  paper  establishes  lower  and 
upper  botinds  on  various  measures  of  Interest  provided  the  demand  rate 
(y)  exceeds  the  production  rate  (c=l) • The  model  presented  in  this 
section  models  the  system  in  continuous  time  on  a continuous  state  space; 
using  results  from  queueing  theory,  analytical  results  are  found  for 
the  evaluative  measures  provided  the  demand  rate  is  less  than  the 
production  rate  (y  < 1) . 


5.  DISCUSSION 

This  paper  has  considered  three  distinct  models  for  studying  inven- 
tory systems  for  perishable  items.  The  key  similarity  across  the  three 
models  is  that  for  each  model  assun^tions  were  made  so  that  the  system 
could  be  characterized  by  a one-dimension  stochastic  process.  For  each 
of  the  models,  analytical  expressions  have  been  found  for  four  performance 
measures:  expected  outdates,  expected  shortages,  average  age  of  an  issued 
iinit,  and  expected  inventory  level.  It  is  felt  that  these  expressions 
can  be  used  for  system  evaluation  and  design  decisions  for  systems  for 
which  the  assumptions  of  the  respective  model  are  reasonable. 
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